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Abstract We consider diffusion-like paths that are ex¬ 
plored by a particle moving via a conservative force while 
being in thermal equilibrium with its surroundings. To 
probe rare transitions, we use the Onsager-Machlup (OM) 
functional as a path probability distribution function for 
double-ended paths that are constrained to start and 
stop at predesignated points after a fixed time. We ex¬ 
plore the continuous-time limit where the OM functional 
has been commonly regularized by using the Ito-Girsanov 
change of measure. When used as a path measure, the 
Ito-Girsanov expression generates an ensemble of double- 
ended paths that are unphysical. We expose the underly¬ 
ing reasons why this continuous-time limit does not, and 
cannot, generate a thermodynamic ensemble of paths. 
Furthermore, we show that the concept of the Most Prob¬ 
able Path and the Thermodynamic action are incompat¬ 
ible with such measures for discrete or continuous time 
diffusion processes. 

Keywords Brownian Dynamics • Pinned Diffusion • 
Thermodynamic Fluctuations 
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1 Introduction 


Much of today’s work in the study of diffusion processes 
is grounded in an expression for the probability of a suc¬ 
cession of states of a spontaneously fluctuating thermo¬ 
dynamic system. This expression, famously reported in 
the 1953 article of Onsager and Machlup 30 , has become 
known as the Onsager-Machlup (OM) functional and has 
so far withstood the test of time. One commentary even 
stated that the results are “incapable of improvement ei¬ 
ther in form or in their mode of derivation” 28 . Here, as 
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in the original article, we represent the fluctuations by 
white noise, spatially and temporally uncorrelated, and 
whose amplitude is given by the fluctuation-dissipation 
theorem. The underlying motion is then expressed in 
terms of Brownian dynamics. 

For the purposes of this paper, we divide rare events 
into two groups: those consistent with thermodynamics 
and those that are not. As an example of the latter, 
what we call extremely rare events, consider the very 
small probability that all the molecules in a room might 
migrate to one corner. Even though this is allowed by 
statistical mechanics, if we observed such an event, we 
would have witnessed a violation of thermodynamics. We 
are not, nor were Onsager and Machlup 30 , interested 


in these extremely rare events. We are interested in rare, 
events that are consistent with thermodynamics; ones 
that are driven by the fluctuations inherent in a thermo¬ 
dynamic system. 

In the 1970s, the original OM functional was ex¬ 
tended to the continuous-time limit using a Ito’s lemma, 
the Girsanov theorem, and the Radon-Nikodyn deriva¬ 
tive. The probability measure, which we will refer to 
as the Ito-Girsanov measure, was written for this limit, 
and has generated considerable attention 


18-21 23 24 35 36 40 41 . This limiting procedure has 


been accepted as a method that can be used to look at 
rare events by constructing probabilities for pinned dif¬ 
fusion paths (5j[l4j[T5j[22][26|[32j[39] . An extension of this 
limiting procedure, proposed by Graham 15 , and later 


Eyink 10 , stressed the generalization of a “least-action” 
principle to describe particle motions, which then leads 
to the notion of the Thermodynamic action. 

In this paper, we show that this Ito-Girsanov mea¬ 
sure, the continuous-time limit of the OM function, is not 
the appropriate probability measure for exploring ther¬ 
modynamically driven transitions. Furthermore, we show 
that the probabilities that arise from such “thermody¬ 
namic actions” are independent of the details of the par¬ 
ticle motions as the probabilities are a result of the noise 
that originates in the thermal bath. Thus there is no 
“action” to minimize. We explain how the Ito-Girsanov 
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functional is unpliysical in that it is inconsistent with 
thermodynamics, and our direct sampling confirms this. 
In addition, for thermodynamic significant events, the 
path-probability measure is flat, that is, it attains the 
same value for any event that is consistent with ther¬ 
modynamics. The point that is missed in the long list 
of works, cited above, is that the Ito-Girsanov expres¬ 
sion only provides the change of measure relative to the 
free Brownian Bridge. Nowhere has it been proven that 
it can be used to give the correct relative probability of 
any two paths. On the contrary, here we show that the 
Ito-Girsanov measure is not a path-probability measure, 
it is only an indicator of how much a free-Brownian so¬ 
lution differs from a solution for a nonzero force. 


2 Brownian Dynamics and the Continuous-Time 
Path Probability 


Throughout this paper, we will consider a particle in con¬ 
tact with a heat reservoir at a temperature e. It is mov¬ 
ing under the influence of a potential V(a;) with the force 
being F(x) = —V'(x). Note that although the equations 
are written for the one-dimensional case for clarity, the 
formalism can easily be extended to higher dimensions 
and for a collection of particles. 

The equation of motion for Brownian dynamics is 
given by the Stochastic Differential Equation (SDE): 

dx = F(x) dt + y/2~edW t (1) 


where dW t is the standard Wiener process that repre¬ 
sents the (uncorrelated) Gaussian noise. Using a discrete 
time step, At, one typically uses the Euler-Maruyama al¬ 
gorithm [25] as an approximate method for propagating 
the position as a function of time. In particular, 

Xi + 1 = Xi + F(xi) At + VYeAt & 


where is a Gaussian random variate with mean zero 
and unit variance. Successive application (N times) of 
this equation produces a sequence of positions {a;*} which 
is called a path. Onsager and Machlup 30 used the un¬ 


derlying thermal fluctuations to write the Gaussian path 
probability, P p oc TUexp(—), in terms of the path 
variables themselves, namely, 


, _ _ At ^ 1 

~ lnP r- C+ U^2 

i=1 


%i +1 

At 


- F (Xi )) 


( 2 ) 


where C is a constant which is unimportant for this pa¬ 
per. This equation defines what is sometimes called the 
thermodynamic action and other times, the OM func¬ 
tional. In the continuous-time limit, using Ito calculus 
and the Girsanov theorem, the Radon-Nikodym deriva¬ 
tive is used to express the change in the measure 
We define the Ito-Girsanov measure to be 
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exp (- (y(x t ) - V(x 0 ) +J dt G ( x t )) 


( 3 ) 


where T is the duration of the path, P p is the continuous¬ 
time limit of P p , Q p is the measure associated with free 
Brownian motion, and the function G(x) is defined as 
G(x) = \F(x) • F(x) - eV"(a). 


3 Away from the Continuous-Time limit 


One of the uses of the OM functional is to incorporate 
it into a scheme to sample paths that are constrained at 
both ends. The aim is to efficiently generate an ensemble 
of paths that include a transition over an energy barrier. 
When the barrier is large compared to the typical ther¬ 
mal energy, the transition is a rare event. We wish to 
explore barrier hopping that is consistent with thermo¬ 
dynamics, where the driving noise reflects the fluctuating 
random effects that originate in a thermal bath and thus 
are independent of the particle’s position. 

The quandary is that for a very simple one-dimensional 
example, the generated paths quickly become unphysi¬ 
cal when using the Ito-Girsanov form, P p , for the path 
probability. Long paths, generated with small time steps, 
are expected to be consistent with equilibrium thermo¬ 
dynamics. First we address the question of what happens 
away from the continuous-time limit and show that it is 
possible to use a method which generates a collection 
of paths that are consistent with the Boltzmann distri¬ 
bution. Then we show that when we use Ito-Girsanov 
expression, given by Eq. [3j we obtain results that are 
unphysical. 

We have constructed a one-dimensional potential to 
highlight one of the key elements in this problem, namely, 
the role of entropy. We explore whether entropic effects 
are properly included by forming a potential having two 
degenerate wells with different widths. In particular, con¬ 
sider the potential (plotted in Fig. [I]) 

Y(x) = 2 -26 (8 — 5x) 8 (2 + 5x) 2 (4) 


which has two degenerate wells with a barrier of unity 
at the origin. A narrow (quadratic) well is on the left 
and the wide well is on the right. By making the wells 
degenerate, we eliminate the large (exponential) depen¬ 
dence due to energy differences. And we accentuate the 
entropic effects by making one well much wider than the 
other. 

To understand how the time discretization affects the 
numerical solution to the SDE we used both Eq. [2] and 
the following form of the path probability P p , 


-inn» p = c+K^ 

i=1 




At 


-F(xi) 


-Ji 


( 5 ) 


Xi is the midpoint 
is related to the 


where Ji = ^fln^l — ^ F'(xt 

of Xi and Xi+\. This expression, P p , 

Stratonovich representation of the diffusion process (see 
Elber and Shalloway [§]). As noted by Van Kampen 38 
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Fig. 1 A plot of the double-well potential used in this pa¬ 
per (see Eq. [4|. A narrow (quadratic) well, on the left, is 
separated from the broad well on the right by an unit energy 
barrier. 

for the problem described by Eq. |T] the limits of both 
the Ito and Stratonovicli discretizations are equivalent. 
In the continuous-time limit Jj => eF'(xi) which pro¬ 
vides the Laplacian term in the definition of G, giving 
the same limit as seen in the Ito-Girsanov expression in 

Eq-i 

In sampling P p and P^, we generated a sequence of 
paths at a temperature e = 0.25, constrained to start in 
the narrow well and to end in the broad well and used a 
HMC implementation using the “Implicit Algorithm” in 
Beskos et al. [I]. 

We use the Heaviside function 0 to define the func¬ 
tion B(s) to be the fraction of the path that is contained 
in the broad well, namely, 

B (s) = ^ dt 0{x[ s) ) « e ( x i S) ) 

where the sampling index is denoted as s, and the cor¬ 
responding path is { }. 

As can be seen in Fig. [5] using either of the path prob¬ 
abilities given by a discrete sum, we find an ensemble of 
paths, which are concentrated in the wide well. These 
paths are constrained to undergo at least one transi¬ 
tion; we find only one transition in the majority of paths. 
We also find paths which undergo multiple transitions. 
Using either expression, on average, the percent of the 
time spent in the broad well is approximately 80% which 
is much closer to the equilibrium value of 90%. This is 
the type of behavior that one expects for paths consis¬ 
tent with the Boltzmann distribution. Away from the 
continuous-time limit, we found a behavior that is close 
to the physical result. In the next section, we show that 
these are in stark contrast to calculations that sample the 
path probability expression given by the Ito-Girsanov ex¬ 
pression (Eq. [3]). 


4 Unphysical Results 

Now we turn the results found when sampling from the 
path probability given by the Ito-Girsanov expression, 
P p (Eq. [ 3 ]). Results generated using the Ito-Girsanov 
expression has been questioned before. For example in 
Adib’s paper l], the origins of the poor performance 
were not understood. The potential, described above and 
used in this work, was designed to highlight entropic ef¬ 
fect on transition paths. With degenerate wells, entropy 
enhances the time the particle will spent in the wide 
well. However, from the form of P p , the curvature of the 
narrow basin would seem to enhance the probability of 
paths that spend a long time there, see the early pa¬ 
per of Weiss and Haffner [40]. Indeed the Most Probable 
Path [31] (MPP) would be highly concentrated in the 
narrow well. 

We used the Ito-Girsanov formula with the method 
developed by Beskos et al 3 , referred to here as path- 
space Hybrid Monte-Carlo (psHMC), using a tempera¬ 
ture e = 0.25, to generate a sequence of paths constrained 
to start in the narrow well and to end in the broad well. 
In Fig. [2j we plot B(s ) using the described procedure 
(red curve), showing that the paths quickly become un¬ 
physical in that the particle spends the vast majority of 
the time in the left, thin well. Such paths are inconsistent 
with the equilibrium thermodynamical distribution. Al¬ 
though only one sampling set is presented here, we note 
that this effect is robust in that similar results are ob¬ 
tained for a range of parameters. 



Fig. 2 The fraction of the path in the broad well is shown as 
a function of an arbitrary sampling index, s. These are results 
for three calculation. We started with the same input path 
where B( 0) « 0.6, and had a path length of T = N At = 150, 
and a time step along the path of At = 0.005. The lower 
(red) curve corresponds to the results using the Ito-Girsanov 
form, while the upper (green) curve are the results of using 
the functional (Eq. [2| with the upper (blue) curve are the 
results of using the functional (Eq.pJ based on a Stratonovich 
discretization. The (orange) dotted line corresponds to the 
equilibrium value. 
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Fig. 3 Results for the calculation using the RWM algorithm 
starting from a randomly generated Brownian Bridge. As a 
function of iteration number, a plot of the fraction of the 
path which resides in the Broad Well (black line) and the 
equilibrium value (orange dotted line). 


5 Random Walk Metropolis 


We now turn to a simpler sampling method, the Random 
Walk Metropolis (RWM) algorithm, to sample P p , which 
proceeds as follows: First, create a Brownian path that 
starts and ends at the origin. A new Brownian Bridge 
is generated and then combined with the current path 
(with the sum of the squares of the mixing coefficients be¬ 
ing unity). A proposed path is generated from the Brow¬ 
nian Bridge by shifting the starting point to be x~ = 0.4 
and the ending point to be x + = 1.6. The acceptance or 
rejection of this new path is based on P p (Eq. 131) . This 
procedure was repeated thousands of times ana gener¬ 
ated the results displayed in Fig. [3] As shown in this 
figure, the starting Brownian Bridge evolves to a path 
that is concentrated in the broad well. The fraction of 
the time that the path spends in the positive half plane 
begins at a value below 0.40 and settles to a value larger 
than 0.95. This result differs from the sampling of the 
Ito-Girsanov expression shown above. However, if one 
starts with a path that is highly concentrated (94%) in 
the narrow well instead of a pinned Brownian path, one 
finds that the path does not evolve significantly. Look¬ 
ing at Eq. [3j one finds that the integral of G over the 
initial path is a large negative number. The mixing pa¬ 
rameter in this case had to be reduced by an order of 
magnitude otherwise any proposed path was rejected by 
the Metropolis step. Clearly, the integral of G derives its 
large negative value from the path being concentrated 
in the narrow well. Looking at the sine transform of the 
path, one would see that the low frequency part of the 
spectrum looks more like that of an Ornstein-Uhlenbeck 
(OU) process 37 rather than that of free Brownian mo¬ 
tion. Thus the proposed moves generated in the RWM 
method are ill-suited for probing the fluctuations around 
this sharp probability maximum. 


6 Numerical Aspects 


To recap the results presented above: we found that by 
sampling pinned paths we were able to generate different 
results. On the one hand, we generated an ensemble of 
paths that were inconsistent with the Boltzmann distri¬ 
bution by using the Ito-Girsanov expression and 1) using 
psHMC algorithm or 2) with the RWM starting from a 
very particular initial path. On the other hand, we gen¬ 
erated paths that concentrate in the wide well, using a 
form of the OM functional with an HMC algorithm 3) 
with an Ito discretization or 4) with a Stratonovich dis¬ 
cretization, away from the continuum limit, and 5) using 
the Ito-Girsanov expression with the RWM starting from 
a Brownian Bridge. 

The key to understanding numerics behind these re¬ 
sults is to examine the “cross term” in Eq. [2] This term 
can be written as 


Xj+i - Xj 

At 


F(xi) 


5±L^.(f(*, +1 )-fW) 

5±a_±l.(F(.T i+1 ) + F(x,)) 


The last term above is an approximation to the change- 
in-potential term in Eq. [3] (as it telescopes under sum¬ 
mation). The previous term is approximately 


%i-\-1 %i 

2 ~At 


(F(x i+1 ) - F(xifj 


Ax 2 
2 ~At 


F'(x) 


( 6 ) 


where x is the midpoint between x t and Xj r+ \. A further 
simplification can be done, but only if the time increment 
is small and the noise and positions are uncorrelated. 
For sufficiently small At , one can average over several 
time steps. Doing so, we arrive at an expression that 
approaches the Ito-Girsanov expression (Eq. [3]) , namely, 


— e Y"(x). (7) 

Remember that the noise, as embodied in Z\a; 2 , must 
be uncorrelated to the positions, as required in the un¬ 
derlying SDE (Eq. [TJ) when probing thermodynamically 
consistent motions. It is important to note that this re¬ 
quirement is satisfied (by construction) when the RWM 
method is used, but only if the correlation is absent in 
the initial path. This means that the some of the symp¬ 
toms of the unphysical probability measure will remain 
hidden as the RWM method lacks the flexibility to fully 
sample the measure. 

In the psHMC algorithm the implied noise history 
and the particle position are intertwined as one searches 
for the paths that have the larger probabilities as indi¬ 
cated by the Ito-Girsanov expression. During the search, 
correlations are introduced into the low end of the fre¬ 
quency spectrum while the path evolves under the (de¬ 
terministic) Hamiltonian flow. High frequencies are cor¬ 
rectly handled as paths generated by both methods have 
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quadratic variations that are close to ideal: J2 Ax 2 ~ 

2 e T. The implication of such low-frequency correlation 
is that the thermal bath is no longer independent of the 
system and thus the fluctuations are no longer thermo¬ 
dynamic, and entropy is generated [9 . 

The resulting paths (not shown here) using the psHMC 
method look like a noisy version of the so-called most 
probable path 31 , MPP. Recall that the MPPs were 


generated by finding the probability maximum implied 
by the Ito-Girsanov expression. Thus we conclude that 
the unphysical nature of the MPPs found previously 311 
is a direct result of using this unphysical measure. RWM 
produces paths that are only mixture of Brownian paths 
which evidently lack the flexibility to probe the regions 
that the measure indicates are of high probability. 

Understanding this cross term (Eq. [6]) helps to ex¬ 
plain why we see less of a problem with calculations for 
an OU process. When F' is a constant, the cross term is 
averaged correctly if the grid spacing is fine enough so 
that the Quadratic Variation is close to ideal. 

We note that our findings are not touched by the 
Ito vs. Stratonovich controversy. In the problem consid¬ 
ered here, the noise and the “diffusion coefficient” are 
spatially homogeneous. As pointed out by Van Kam- 
both discretization schemes lead to the same 
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correct limit. Looking at Eq. [5j it is clear that in the 
limit as At —> 0, we recover the Ito-Girsanov functional, 
Eq. [3j if we regularize the expression. It is the regu¬ 
larization step that hides an assumption. When using 
the Radon-Nikodym derivative, the Ito-Girsanov mea¬ 
sure cannot be the probability measure because it has 
meaning only when using free Brownian Bridges and not 
for a general path. 


7 Equivalence of Path Probability 

We now turn to a more general discussion of the idea 
of an MPP in the OM formalism. We first examine the 
relative probability of paths. As in any continuous dis¬ 
tribution, while the probability of any one path is zero, 
the relative probability of two paths is well-defined. In 
the case of free Brownian motion, for sufficiently long 
paths of the same duration, the relative probability of 
any two paths is unity. But this is also true when con¬ 
sidering Brownian dynamics with any force. In all cases, 
the path probability is governed by the Gaussian noise, 
P p oc lit exp(— |£?), where the set {^} is independent 
of the force and thus independent of the position of the 
particle. The same argument applies to any finite repre¬ 
sentation of the path, and thus holds as the size of the 
time increment becomes infinitesimally small. 

To illustrate the above point, consider the following 
“thought” experiment. For a conservative force consider 
sampling the Boltzmann distribution using the Brown¬ 
ian dynamics as expressed by Eq. |T] Take the starting 
point, xq, to be arbitrary and integrate the SDE over 



Fig. 4 Distribution of endpoints of the set of trajectories 
with equal probability. Here the endpoints of 472,640 differ¬ 
ent trajectories are shown. The blue line is the normalized 
Boltzmann distribution, exp(—Y/e)/Z, for the temperature 
e = 0.25, and potential given by Eq. [4] 


a fixed time, T, that is long compared to any barrier 
hopping time. Using a non-zero but small time step, At , 
one uses N r = T/At Gaussian random variates. Keeping 
the same set of random numbers, but simply scrambling 
the order, redo the integration. This provides N r \ paths 
each with identical probabilities. For large enough T and 
small enough At, the set of endpoints, x(T) should be 
distributed in a manner that is close to Boltzmann. Here 
it is important to recognize that it is not the path proba¬ 
bility that creates the distribution, as all paths are iden¬ 
tically probable. Rather it is the path density that drives 
the correct distribution of the endpoints. We have par¬ 
tially preformed this experiment (calculation) for the po¬ 
tential defined in Eq. |4j We chose N r = 2 million random 
numbers, Gaussian distributed with mean zero and vari¬ 
ance one. Out of the possible N r \ paths we followed over 
400,000 of them using the Euler-Maruyama algorithm, 
with At = 0.0005. The endpoints are then collected and 
the histogram is displayed in figure [4j Evidently the dis¬ 
tribution looks to be very close to the Boltzmann distri¬ 
bution in spite of the fact that each path has identically 
the same probability. Similar results are expected when 
one generates independent sets of random numbers in¬ 
stead of using the permutations of the original set. 

The noise is a consequence of the random fluctua¬ 
tions of the thermal reservoir, and in the SDE given in 
Eq.0 the noise is not correlated with the position of the 
particle. These considerations point to the inappropri¬ 
ateness of using the Ito-Girsanov change of measure as 
the path probability functional for generating thermody¬ 
namically consistent paths. By inspecting Eq. [3j we see 
that it gives different probabilities for different paths: 
some paths are more probable than others. This is ac¬ 
complished by correlating the noise with the positions 
through the function G. It is clearly incorrect to use the 
Ito-Girsanov expression in this manner. And in particu- 
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lar, this is the immediate origin of the numerical results 
displayed above. 


8 The Meaning of The Change of Measure 


The question is how do we understand our results in light 
of the above information? To accomplish this, we look at 
the ramifications for an Ornstein-Uhlenbeck (OU) pro- 
with the force being Fqu = —"/x. The Ito- 
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Girsanov change of measure indicates that the probabil¬ 
ity distributions of two diffusions that differ only in their 
drift term are mutually absolutely continuous. Exact so¬ 
lutions are known for both the free Brownian motion and 
the OU process. For the OU process, the frequency {v) 
spectrum is finite at the origin with al/v tail. The fre¬ 
quency spectrum of the free particle motion falls off as 
1/v after diverging at the origin. For a long enough path 
length, every realization of the Brownian Bridge has this 
same spectrum. Since linear combinations of such realiza¬ 
tions are still (free) Brownian Bridges, they too have the 
same spectrum. Thus each Brownian Bridge and every 
linear combination of Brownian Bridges have a frequency 
spectrum that differs from that of an OU process; they 
are not solutions to the OU SDE. This holds true for 
other forces; in general, in the limit as the path length 
becomes infinite, free Brownian Bridges are not solutions 
to the SDE (Eq. |T|) with nonzero drift. This then explains 
why the RWM sampling method is increasingly sluggish 
as the path length increases. 

With this in mind, we now examine the Kullback- 
Leibler divergence Dkl(P\\Q) between the measures for 
two diffusions P and Q that differ in their drift terms. 
If we take Q to represent the measure of free Brown¬ 
ian Bridges then the expectation of the Ito-Girsanov ex¬ 
pression can be identified as the Kullback-Leibler diver¬ 
gence since both measures have the same “normaliza¬ 
tion” terms as inferred from the discussion in the pre¬ 
vious section. This then explains why the Ito-Girsanov 
change of measure should be viewed as an indicator of 
the differences between the solutions of two SDEs that 
differ only in their drift term and is not related to the 
probability of paths of one or the other diffusions. 


9 Conclusion 


The Onsager-Machlup functional is based on Brownian 
dynamics and provides a way of understanding the dou¬ 
ble ended path sampling problem 15 . It had been ac¬ 


cepted that in the continuous-time limit, the Onsager- 
Machlup functional could be replaced by the Ito-Girsanov 
change of measure (Eq.|3| as a way of handling the infini¬ 
ties inherent in such a limit. We have shown here in three 
ways that this is inappropriate. First, direct sampling 
gives unphysical results. Second, in the long path limit, 
any linear combination of free Brownian paths cannot be 


a solution to the OU process nor to diffusions with more 
complicated forces. Third, interpreting the Ito-Girsanov 
change of measure as a probability distribution favors 
some paths over others even though paths of the same 
duration must have the same probability. We determined 
that the Ito-Girsanov change of measure is simply an in¬ 
dicator of how inaccurately the free Brownian paths rep¬ 
resent solutions to the diffusion process with a nonzero 
drift. Thus using the Ito-Girsanov change of measure as 
a probability distribution is not thermodynamically cor¬ 
rect. 


The sampling of paths is a complicated numerical 
procedure interlaced with sophisticated mathematics. In 
such problems, it can be difficult to separate concep¬ 
tual errors from numerical artifacts. It was previously 
observed that the MPPs [31] for many potentials seemed 
to be unphysical, which should have been taken as an in¬ 
dication that something was amiss. Now we understand 
why. The Ito-Girsanov expression is not a path probabil¬ 
ity measure and thus the minimizers of the Ito-Girsanov 
expression are not related to physical paths. 

The particular potential, used here, was constructed 
to highlight entropic effects on transition paths. Using 
it, we were able to pinpoint a common misconception. 
We have found similar behaviour with other potentials, 
such as two quadratic wells with differing widths (con¬ 
nected with a third-order spline). Also, in two dimen¬ 
sions, we studied diffusion in the potential V (x, y) = 
(x 2 +y 8 — l) 2 exp(— 2x 2 — 3(y + l) 2 ), which has an almost 
degenerate channel. There we found unphysical results 
for a temperature e = 0.10; these results are consistent 
with the MPP in that the path has a tendency to con¬ 
centrate near the point (0,1) where the channel walls are 
steepest, rather than in the broad wells. But such details 
are unimportant in that they again demonstrate similar 
behavior as found in the one-dimensional case that we 
have reported here. The results for the simple double¬ 
well potential shown in this paper clearly illustrate that 
the Ito-Girsanov expression does not correctly describe 
entropic effects when used as a basis for sampling paths. 

The findings in this article affect a wide body of 
work. As mentioned above, the theory for entropy pro¬ 
duction in nonequilibrium thermodynamics, such as the 
work of Speck et al 34 , should be re-examined. In¬ 


deed it would seem that even an encyclopedia entry 27 


needs to be modified. For probing the folding of pro¬ 
teins, the recent work of Fujisaki, et al 13 suffered from 
using the unphysical form of the Onsager-Machlup func¬ 
tional. It is unclear how to interpret their results. Other 
works |7|11|17|42| have to be reevaluated in light of these 
new results. 


The results reported here also have ramifications for 
sampling general measures in infinite dimensions. The 
samples or paths in infinite dimension must be constructed 
to probe the high probability regions as defined by the 
measure. Algorithms based on Brownian Bridges have a 
very limited flexibility in that the noise history is un- 












7 


correlated by construction. Only in HMC methods, high 
probability regions are explored due to the Hamiltonian 
flow which allows the high probability, highly correlated, 
modes to be explored. The limited sophistication of the 
sampling algorithms employed in the past has hid some 
of the results exposed in this article. 

The path integrals that appear in diffusion look sim¬ 
ilar to those that appear in Quantum Mechanics. The 
natural question is what is impact of the current work on 
Feynman Path Integrals 12 . The short answer is none. 


Firstly the paths in Quantum Mechanics are smooth, 
while Brownian paths are almost nowhere differentiable. 
Secondly diffusion is a local process while solutions to the 
Schrddinger equation are not. The time evolution of the 
wave function is most easily found by using an eigenfunc¬ 
tion expansion; the eigenfunctions are dependent on the 
potential everywhere in space. Diffusion is described by 
a Markovian process and white noise; local in space and 
time. Thirdly, in the path integral formulation of Quan¬ 
tum Mechanics, the path distribution is peaked around 
the classical path, a smooth curve that is defined by the 
classical action. As we discussed above, each physical 
(thermodynamically-consistent) diffusion bridge has the 
same probability. The most important observation in this 
regard is that the field-theoretic methods, such as in the 
review by Smith [33], developed for Feynman integrals 
are not transferable to classical diffusion. 
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